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Discrete breathers (nonlinear localised modes) have been shown to exist in various nonlinear Hamil- 
tonian lattice systems. In the present paper we study the dynamics of classical spins interacting via 
Heisenberg exchange on spatial d-dimensional lattices (with and without the presence of single-ion 
anisotropy). We show that discrete breathers exist for cases when the continuum theory does not 
allow for their presence (easy-axis ferromagnets with anisotropic exchange and easy-plane ferromag- 
nets). We prove the existence of localised excitations using the implicit function theorem and obtain 
necessary conditions for their existence. The most interesting case is the easy-plane one which yields 
excitations with locally tilted magnetisation. There is no continuum analogue for such a solution 
and there exists an energy threshold for it, which we have estimated analytically. We support our 
analytical results with numerical high-precision computations, including also a stability analysis for 
the excitations. 
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I. INTRODUCTION 

The phenomenon of dynamical localisation has been 
a subject of intense theoretical research. It is well 
known that classical Hamiltonian lattices possess peri- 
odic in time and localised in space solutions called dis- 
crete breathers or intrinsic localised modes. A recent ex- 
plosion of interest to discrete breathers has happened due 
to the fact that they may exist in lattice models of inter- 
acting identical 'pavticles. Breathers in continuum models 
(like, for example, the well known sine-Gordon equation) 
exist only due to the high symmetry of the system and 
are therefore structurally unstable. Discrete breathers 
are generic solutions of nonlinear lattice equations. Their 
existence is based on the fact that all possible resonances 
of multiples of their frequency with the bounded band 
of small amplitude plane waves (BSAPW) above the 
classical ground state can be avoided. So far discrete 
breathers have beeM proven to be generic solutions in 
both HamiltonianaH and dissipativea systems. Several 
cases of experimental observation of discrete breathers 
have been reported [in Josephson junction arrays, ar- 
rays of|-,weakly coupled waveguidesO, low dimensional 
crystala3, and biological systems (myoglobin)Ll| . 

Due to its spatial periodicity, lattices of interacting 
spins are ideal systems to observe discrete breathers as 
well. Here we will concentrate on large spins which may 
be described classically. Nonlinear waves in magnetic sys- 
tems hai/e-been extensively studied during the last three 
decadeaj'EHl. The results of these studies provide a lot of 
information about the properties of solitary waves (par- 
ticularly, breathers) in magnets, since it is possible in 
many cases to obtain explicit solutions to them. How- 



ever, neglecting discreteness effects may lead to loosing 
important features of the nonlinear wave dynamics. For 
instance, since only high symmetry continuous systems 
possess breather solutions, the area of the potentially in- 
teresting models is artificially reduced. Another draw- 
back of continuous systems is that the consideration of 
nontopological localised excitations is typically restricted 
to space dimension one. 

In the last decade, a number of papers has appeared 
where localised modei s i n .magnets were treated as essen- 
tially discrete objectstilll3 (also an attempt of experimen- 
tal observation-af discrete breathers in antiferromagnets 
has been madeO). However, no rigorous existence proofs 
have been given, and only the most simplest cases (from 
the point of view of symmetries) have been considered. 
Preserving these symmetries one can continue those so- 
lutions to the space continuous limit. 

The aim of this work is to present breather excitations 
for spin lattices in cases where the symmetries will not 
allow for a similar mode construction in space continu- 
ous cases. Also we will not restrict the consideration to 
one-dimensional systems. Below we present a rigorous 
existence proof for discrete breathers in magnetic sys- 
tems using the anti-continuum limit. With the help of 
this proof and of a Newton iteration method, we show 
the existence of discrete breathers in ferromagnetic lat- 
tices with anisotropic exchange interaction. We also con- 
sider an easy-plane ferromagnet and find a new type of 
discrete breathers, with several spins precessing around 
the hard (single ion anisotropy) axis, while all the oth- 
ers precess around an axis which lies in the easy plane. 
Note that only the simplest case of monochromatic in 
time breathers has been investigated in most of previ- 



ous papers. Such a situation simplifies considerably the 
treatment of the system, and important families of solu- 
tions can be lost. Our studies do not depend in any way 
on the number of higher harmonics in the time evolution 
of a breather. 

This paper is organised as follows. The next section 
presents the model Hamiltonian and the equations of mo- 
tion. In Sec. Ill, we consider an easy-axis ferromagnet, 
discuss the implementation of the anti-continuum limit, 
and give a rigorous existence proof for discrete breathers. 
In Sec. IV, we study an easy-plane ferromagnet. Then 
Sec. V presents a study of a two-dimensional lattice with 
easy-plane anisotropy. Discussions and conclusions arc 
given in Sec. VI. 



anisotropy {Jx,Jy « Jz), or by introducing an on-site 
anisotropy term D > 0. Before we study breather so- 
lutions of Eqs. (||), let us study the dispersion laws for 
linear spin waves. 



A. Dispersion law^s 

First, we consider the easiest case, a ferromagnetic 
chain without ion anisotropy (Z? = 0), but with a strong 
exchange anisotropy: < Jx, Jy << Jz- In this case, 
the ground state is S'f^ = ±1,5',^ = 3^ = 0. Linearis- 
ing the equations of motion around one of these ground 
states, e.g., 5"^ = Sx sin {qn — ujt), S^ = Sy cos {qn — ujt), 
S^ — const « 1, we obtain 



II. HAMILTONIAN AND EQUATIONS OF 
MOTION 

We consider a lattice of classical spins described by the 
Hamiltonian with Heisenberg XYZ exchange interaction 
and single-ion anisotropy 



^lil) = (Jz - Jx COS q){Jz - Jy cosq) 



(4) 



H 



n^^n' a = {x,y,z) n 



Here S^.Sl^, S^ are the nth spin components (n labels a 
lattice site) which satisfy the normalisation condition 



qx2 I 01/2 I qz2 _ n2 



(2) 



For simplicity, the total spin magnitude can be nor- 
malised to unity: S — 1. The constants Jx,Jy,Jz > 
are the exchange integrals and D is the on-site anisotropy 
constant. 

The equations of motion for the spin components in the 
one-dimensional spin chain with nearest-neighbour inter- 
actions are the well known Landau-Lifshitz equations: 



^Ds^s:, , 



(3) 



^n — 2 W^'^n \^n~l + ^n+l) Jy^n i'-'n-l + '-'ri+ljj • 

Generalisation to higher lattice dimensions is straightfor- 
ward. 



III. EASY-AXIS FERROMAGNET 

First, we consider spin lattices with a ground state 
that corresponds to all spins directed along a given 
axis (we assume this axis to be the Z-axis). This can 
be achieved either by introducing a strong exchange 



This dispersion law is shown in Fig. 1, with the edges of 
the linear band uj^ (?) being given by 



\'-' z ^x}\'Jz '^yji 



{Jx+Jz){Jy + Jz) ■ (5) 



C5^ 




nil 



FIG. 1. Dispersion law for the ferromagnetic chain with 
strong exchange anisotropy. 

Let us explain how we can use the information about 
dispersion relations in order to formulate expectations 
about existence or nonexistence of breathers. It has been 
shown inliJ that breather solutions of the full nonlinear 
equations bifurcate from certain nonlinear plane wave so- 
lutions. These specific plane wave solutions are time pe- 
riodic and reduce to band edge plane waves (BEPW) in 
the limit of small amplitudes, i.e., to linear waves with 
^L{q_BEPw) being an extremum of u!L{q)- A necessary 
prerequisite of the existence of breathers is that the fre- 
quency of these BEP W's is detuning away from the linear 
band (^^(g) with increasing amplitude or energy density 
of the wave. 

Let us consider an excitation which corresponds to the 
value of the wave number q = Q from above. This is a 
space-homogeneous excitation S^ — S" and the Landau- 
Lifshitz equations yield 



5" - {Jy - J.)sys' , 
sy = (J, - j.,)s^s' , 



(6) 



Jy we 



As can be seen from these equations, even for X 

have softening in the dispersion law at g = 0. Then S^ — 

0, and for 5"^ = A cos uit the spin precession frequency 



OJ 



(l - A^){J, ~ J^f < {J.^J^Y 



(7) 



is below the lower edge of the linear band, which indicates 
an occurrence of discrete breathers in the band gap. Note 
that for the completely isotropic model J^ = Jy = Jz = 
J, there is no gap and consequently no breather solutions 
are to be expected. 

It is easy to check that a similar analysis of the upper 
band edge yields lowering of the BEPW frequency with 
increasing amplitude, i.e. its frequency is attracted by 
ujL{q) instead of being repelled. Consequently we do not 
expect breathers to bifurcate from the upper band edge. 

For D j^ there are no qualitative changes. The 
ground state of the chain remains to be the same, whereas 
the gap in the dispersion law widens: 



^£(9) == (Jz - Jx COS q){Jz - Jy cos q) 



AD\Jz-'^^-^cosq 



AD^ 



(8) 



with 



^0 = {Jz - Jx){Jz - Jy) +AD[Jz- •LlLLl \ + 4i52 



(9) 



LOI = (J^ + J,){jy + J,) +4.D[Jz + 



J X I J\t 



AD' 



Note that the band gap exists even in the case of isotropic 
exchange (all Jq's are equal). Therefore, easy axis 
anisotropy increases the chances of breather existence. 



In this case the Z-component of each spin is conserved. 
The solution of (^) reduces to the precession of decoupled 
spins around the Z axis with frequencies which depend on 
the realised value of the Z-components of nearest neigh- 
bour spins (due to nonzero Jz): 



Sn + iSy = Ane'^'^"'+'^-^ , ^,t = 0, 



(10) 



where the precession frequency of the nth spin is given 
by 

\^n\ = y (5,l_i + 5,1+1) + 2DS^, Al = \- Sf . (11) 



1. Vanishing single-ion anisotropy (D — 0) 

The initial choice of one precessing spin and all the 
others being at rest 



1,1,1,^0,1,1,1,...) 



(12) 



with ^0 < 1 cannot be used to generate breathers because 
the frequency of this solution oj = Jz resonates with the 
linear band [more precisely, with the lower edge of the 
linear band loq, see Eq. (|^)]. Therefore it cannot be 
continued to the region of non-zero Jx and Jy. A way 
out is simply to excite three neighbouring spins: 



^n — (..., 1, 1, 1, "S*!, Sq, Si, 1, 1, 1, . 



(13) 



Here < S'o < S*! < 1 and, since the precession frequency 
of all the three central spins must be the same, Sq = 
2Si — 1. In this case, the precession frequency u> = JzSi, 
allowing for the absence of resonances with the linear 
spectrum frequency loq = Jz, must satisfy the condition 



2. Additional single-ion anisotropy (D > 0) 



B. Implementation of the anti-continuum limit 

Now, following MacKay and AubryEJ, we apply the 
anti-continuum (AC) limit to our system. The principle 
of the AC limit consists in decoupling the lattice sites 
and exciting only one or a small number of them, keeping 
all the other in the ground state. Then, upon switching 
on the interaction, the persistence of the localised solu- 
tion is shown. As a prerequisite for the successful exis- 
tence proof and continuation of the breather solution, the 
initial 'decoupled' periodic orbit must be anharmonidi3, 
and the breather frequency and all its multiples should 
not resonate with the linear magnon band. In the case 
of strongly anisotropic exchange {Jx,Jy « Jz), the 
particular case of the AC limit means simply setting 

Jx Jy U. 



In this case we may use the AC limit with the ansatz 
(|l2|). The initial distribution of Z-spin components are 
chosen as follows 



5^ = (...,1,1,1,^0,1,1,1,...) 



(14) 



Here one central spin is precessing with frequency lu — 
Jz + 2DS^ . All other spins are at rest. Small devia- 
tions from their equilibrimn states yield precession with 
frequencies 

wi=wo + ^(S'o-l) , ujo = Jz{l + 2D) , (15) 

distributed along the lattice in the following way 

LOn = (...,Wo,Wo,Wi,W,tJi,CJo,tJO,...) ■ (16) 



Discrete breathers can be continued from the AC hmit if 
the following non-resonance conditions are satisfied: 



kuj ^ ciJo, kiij 7^ tJi, k (z Z 



(17) 



Taking into account that Sq — {uj — Jz)/2D and substi- 
tuting it into the non- resonance condition w ^ wi, we 
get 



{W \4D 2 AD 



(18) 



Note that for fc = 1 the resonance will occur for any 
breather frequency if Jz = 4D. For this set of param- 
eters, a breather continuation from the AC limit is not 
possible for any frequency. 

For this particular case we try another ansatz, namely 
the even-parity case: 

^n — (•■•7 Ij Ij 1) ^0, So, 1, 1, 1, ■•■) , 
W„ = {...,U!o,UJo,UJl,UJ,UJ,UJl,UJo,UJo, :■) (19) 



(20) 



with ojQ and lji being the same as in Eq. (|15|), and 
^ = y (1 + ^o) + 2DSo . 



Then, using the non-resonant condition ( p^ ) which is 
valid for this ansatz as well, we obtain 

kuj i^ ijj\ = ; —UJ . (21) 

^ 2{Jz+AD) Jz+4.D "- ' 

It follows from this expression that the even-parity AC 
limit allows for a continuation of the breather solution 
for all values of Jz and D. 



3. Isotropic exchange Jx 



Jy J z 



J 



Here the ansatz (y^ can be used and the frequencies 
in the AC limit will be distributed as a;„ — 2DS^. The 
eigenfrequencies of the non-excited spins do not depend 
on the values of the adjacent spins and equal ujq- Thus, 
the only non-resonance condition to be fulfilled is koj ^ 

Wo- 



C. Existence proof for magnetic breathers 

Here we present a rigorous proof of the existence of 
discrete breathers for the particular case of strongly 
anisotropic exchange and D > 0. 

Theorem 1 If a periodic orbit of the Hamiltonian ^ 
with a frequency to is non-'o^sonant (kuj ^ Wo,i, k G Z 
and uJL{q)) and anharmoniS^, then the periodic orbit of 



the equations of motion (Hj at a ^ {Jx,Jy} = given 
by the spin precession frequency ^I^ has a locally unique 
continuation as a periodic orbit of the equations fl) with 
the same period T = 27: /lj for a sufficiently small a. 

Proof. Let SLt be the space of bounded infinite se- 
quences z = {zn}nez of triplets z„ = (5,^, S^, S^) of con- 
tinuously differentiable functions of a period T — 2t: /lu 
with symmetry properties: 

s^it) = s:i~t), syjt) = -sy{~t), s^, = 5^(~i) . (22) 

Then the size of oscillations on the nth site will be mea- 
sured by the following norm: 

\zn\ = snv[\Smi \Sl{t)l \K(t)l \SHt)[Ae r}. (23) 

Next, the size of z G SLt is given by 

|z|^sup{|z„|;neZ} (24) 

and therefore SLt is a Banach space. 

Consider now another Banach space SMt of bounded 
infinite sequences w = {wn\nez of triplets w„ = 
( Af,^ , My , M^ ) of continuous functions of a period T with 
the following symmetry properties: 

M:(i) = -M:(-i), Ml{t) = Ml{-t), M^ = -M^i^t) , 

(25) 

and the norms 

\wn\ = sup{|M„-(f)|, |Af„^(t)|, |M^(t) - l|;t e R}, 
|u;| =sup{|ii;„|;neZ}. (26) 

Define the mapping F: SLt — > SMt 

F{z, a) — w 

with a = {Jx,Jy}, z = {S'i^,S'^,S'^}„gz, w = 
{MiJ,M^,M^}„gz, and 

'^n = ■^ \JySn i'-'n-l + '-'ri+lj ~" JzS^ ('-'n-1 + "Jn+ljJ 



2L>j,nS^ D„ , 
1 



(27) 



-'^n — ■^ [JxSn (>S„_i + '-'n+lj Jy^n ('-'n-1 + '-'ri+ljj 

The symmetric solutions of the equations of motion (y) 
are in one-to-one correspondence with zeroes of F, i.e., 
F{z, a) = 0, especially in the case when a = 0. Using 
the implicit function theorem, we prove this solution to 
have a locally unique continuation z{a) for sufficiently 
small a, such that F[{z{a), a] = provided F e C^ and 
its derivative with respect to z, V, is invertible at a = 0. 



To show the invertibility of VF, we linearise the map 
around the periodic orbit at a = 



^M = VFSS 



(28) 



and show that it is invertible simultaneously in the fol- 
lowing three parts of the lattice: the central site n = 0, 
its adjacent sites n — ±1, and the remainder of the lat- 
tice. Invertibility of VF is equivalent to the invertibility 
of the corresponding matrix in (|2S 
(i) For n 7^ 0, ±1 we have 



SM^ 



'SSIloq - 5Sl 



6M^ = -6S^ . 



(29) 



Since the functions SM^ and 5S" are periodic in time, 
we can expand them into Fourier series 



+ OC 



Y^ S...{k)e 



ikut 



(30) 



k— — c 



The time-symmetry requirements (22) ensures S...{k > 0) 
to be either purely real or imaginary, so that 



SM^ii-k) = -5MZ{k) , <5M^(-fc) = 5Ml{k) , 

SM^li-k) = -SM^ik) , 
SS^i-k) = 6S:{k) , SSy{-k) = -SSJiik) , 
<55^(-fc) = 55^ (fc) . 



(31) 



Particularly, SM^{0) = SM^{0) 
result, for fc > we obtain 



6Sy{0) = 0. As a 



SM^{k) = -ujoSSy{k) - ikujSS^ik) 
5Ml{k) = uJaSS^ik) - ikw5Sl{k) , 
5M^(fc) = ^ikojSS^ik) . 



(32) 



These equations appear to be decoupled with respect to 
k and they can be inverted if fc^w^ 7^ loq. The third 
equation can be inverted for k j^ Oiiui ^ 0. For fc = the 
inversion is impossible, but this degeneracy can be lifted 
by imposing the normalisation condition (y). In fact, the 
third equation can be dropped because S^ is defined by 
S^ and 5*^ with accuracy up to a sign. All that we have 
to check is the inequality S*^^ -t- S^'^ — 1 < 0. 

(ii) For n = ±1 we act similarly to the previous case: 



SMf = SSfuji - SSf , 



(33) 



A similar condition for invertibility can be found: k lj ^ 

(iii) Case n = 0. Using the canonical coordinate y and 
momentum p 



we obtain the following Hamilton equations: 

. dH . _ dH 
dp ' dy 



y 



(35) 



We define the pair (w, v) instead of the set {M^ , Mq, Mq ), 
with 



(36) 



dH dH 

u = ^ — y ^ ^ ^ ^p^ ^r ■ 

op oy 

The new function satisfies the following symmetries: 

y{t) = y(-t), p = ujt + h{t), h{t) - -h{-t), 
y{t + T)^y{t) , h{t)^h{t + T) . (37) 

Here h{t) is a periodic function in time. Keeping in mind 
that Sti and SI are fixed by variations of S^f, we obtain 



5u = 2D5y ~Sp^ 2DSy - Sh 

Sv — —Sy . 



(38) 



After considering the corresponding Fourier series with 
respect to time, we find the following inversion condi- 
tions: w 7^ for fc 7^ and Z? 7^ for fc = 0. 

Once the invertibility is shown, by the implicit function 
theorem, the initial solution 2(0) has a locally unique 
continuation z{a) in SLt, and thus the theorem has been 
proved. □ 

This result can be easily extended to lattices in higher 
dimensions, antiferromagnets, and systems with larger 
interaction radius. 



D. A method of computation of discrete breathers 
and linear stability analysis 

For numerical simulations it is convenient to use stere- 
ographic coordinates. The new coordinates incorporate 
the normalisation condition and reduce the problem with 
three unknown real functions per site to the problem with 
one unknown complex function per site: 



e« 



QX 



tsy 



i + s^ 



(39) 



The inverse transform is given by 

ox Sn ' Sit qy _ ^ S» Srt qz 

^ - I- "^' n ■-] I 1^ 12' n 



u 



l + le« 



i 1 + le« 



l + le« 



(40) 



In these new coordinates the Landau-Lifshitz equations 
take the form 



6. 



1 

4i 



[Jx \ fJy) 



Vm-1 SnSn-1 ^ Sn+1 

1 + 16.-1 P 



SriSn+1 



I+I6.+ 1P 



+ (J^ - Jy) ( \ I. ;' + ., I. ?r ) (41) 



S^ = v/r^7cosp, SI = 71^7 sinp, S^^ = 2/, (34) ^J.^n ^^^ |^^_^|2 + ^ ^ |^; 



I + I6.-1P l + |^„^ 



;n-l-l| 



%D^, 



1-I6.P 



l + l?« 



The computation of discrete breathers is done using the 
Newton mapa. 

The hnear stabihty analysis of discrete breathers is per- 
formed by hnearising Eqs. (|l]): ^„(i) = £,n\t) + e„(t) 
around the breather periodic orbit, and solving after- 
wards the eigenvalue problem 



Re e„(T) \ ^ /^ Re e„(0) 
Im e„(T) ) -^ ' I Im e„(0) 



(42) 



If the eigenvalues A of the Floquet matrix Ai are 

found to be located on the unit circle of the complex 
plane, then according to the Floquet theorem, the peri- 
odic orbit is stable, otherwise it is unstable. 



E. Breather solutions in an easy-axis ferromagnet 

Breathers in magnetic lattices with an easy-axis 
anisotropy can be viewed as localised spin excitations 
with the spins precessing around one of the ground states 
of the system (which we have chosen S"^ = 1 at the be- 
ginning of the paper) , so that the effective radius of this 
precession decreases to zero at n ^ ±cxd. The case of an 
isotropic exchange in XY ( J^, — Jy) is the simplest one 
because the S^ component is conserved in the solution, 
and therefore the separation of the time and the space 
variables 



Oj-^ — o^-^ -r "^^n — ^n^ 



(43) 



is possible in the Landau-Lifshitz equation [this can eas- 
ily be seen also from Eq. (pl|)]. The precession ampli- 
tudes An do not depend on time and the Landau-Lifshitz 
equations (|3|) are reduced to a set of algebraic equations 
which can be solved by a simple iteration procedure. This 
is true both in the case of strong exchange anisotropy 
Jx = Jy '^ Jz and also when the exchange is isotropic 
and an on-site anisotropy D > is present. 

The breather existence in the Jx.y — D plane is gov- 
erned by the non-resonance conditions given in Subsec. 
IIIB for one single harmonics k = 1. With the growth 



of J the breather frequency may hit the linear spectrum 
which marks the boundary of breather existence on this 
plane. The nature of the other non-resonance condition 
(O) is different, e.g., we cannot continue a breather solu- 
tion for small J = Jx,y when J^ = 4D, however breathers 
exist for larger values of J. Note, that in the case of 
Jx = Jy discrete breathers have a continuum equivalent 
which is a breather solution of the integrable nonlinear 
Schrodinger equation. The reason is that the XY ex- 
change symmetry allows one to find solutions which are 
monochromatic in time (^3|). As long as the linear band 
provides a gap and the nonlinearity allows for pushing 
the breather frequency into the gap, localised excitations 
may be found regardless of the degree of discreteness of 
the system, which can be characterised by the ratio of 
the gap to the bandwidth of the spin wave spectrum. 



Now we consider a rhombic chain with Jx ^ Jy < Jz- 
Breaking isotropy in the XY plane implies that S^ is not 
conserved in the solution anymore, and according to the 
Landau-Lifshitz equations, it is impossible to represent 
the breather solution in the form (Esl). This implies that 
breathers will have an infinite number of harmonics in 
time, and consequently the spin dynamics is more com- 
plicated. Each spin now draws an 'elliptic' trajectory on 
the unit sphere, elongated towards the larger component 
of Jx or Jy. Fig. I 



^ 




FIG. 2. Dynamics of the central spin with site number 
no = 11 in the chain oi N = 21 spins with Jx = 0.1, Jy = 0.23, 
J, = 1, _D = 0. 



shows such a dynamics of the central spin no = 11 of 
the breather in a chain consisting of iV = 21 spins. The 
breather profile at some instant of time is shown in Fig. 
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FIG. 3. Discrete breather profile at a fixed instant of time 
for the parameters described in Fig. 

Due to the broken symmetry in the XY plane we are 
not able to find breathers in the corresponding contin- 
uum problem. The reason simply is that the linear band 
of a continuum equation may still have a gap but will be 
unbounded from above. Consequently there will be un- 
avoidable resonances of higher harmonics of a breather 
with the linear band causing in general nonexistence of 
the solution itself. Here we have a nontrivial case where 
the discreteness of the lattice provides the necessary sup- 
port for breathers which is missing in the continuum case. 
The computation of breather periodic orbits in this case 
cannot be reduced to solving a system of algebraic equa- 
tions and we have to work in the fiiil phase space by 
using, e.g., a generalised Newton mapEJ. 

Let us briefly discuss the stability of the obtained 
breather solutions. Previous stability studiesllZl have 
shown that the stability depends on the breather parity 
(i.e., its spatial symmetry). Our Floquet analysis of the 
eigenvalues of the stability matrix confirms these find- 
ings. We obtain that the site-centred breathers (contin- 
ued from the one-site breather) are stable in the limit of 
small exchange (see Fig. 0a) 



< 




< 




(b) 



-1 -0.5 



0.5 
Re A 



1.5 



FIG. 4. Eigenvalues {A} of the Floquet matrix for (a) 
site-centred breather and (b) bond-centred breather in the 
easy-axis chain of A'^ = 32 spins with D = 1, J = 0.01, 
u = 1.3. Arrows show the direction of motion of eigenval- 
ues when J is increased. 

while the bond-centred breathers [continued from the 
two-site breathers, see ansatz Eq. (|9|)] are unstable arbi- 
trarily close to the AC limit with the unstable eigenvalue 
being located on the positive half or the real axis outside 
of the unit circle (see Fig. 



IV. EASY-PLANE FERROMAGNET 



In the case of an easy-plane anisotropy we choose D < 
and Jj. — Jy — J^ = J. The ground state of the system, 
without loss of generality, can be assumed to be 



QX 



1, S — b — . 



(44) 



Note that the ground state is degenerate, so that the 
spins can be oriented arbitrarily in the XY plane, but 
they must stay parallel to each other . 



FIG. 5. Schematic representation of a discrete breather 
with one 'out-of-plane' spin in the easy-plane magnet. 

or for two parallel precessing 'out-of-planc' spins, as 
shown in Fig. M. 



A. Linear dispersion law and the anti-continuum 
limit 

Linearising the equations of motion in the vicinity of 
the ground-state (H) we obtain the following dispersion 
law: 



u?{(l) = J^(l - co%qf + 2J\D\{\ - cosq) 
This is an 'acoustic'-type dispersion law with 



(45) 



tj2^w2^0)=0, ul=Lu^{Tr)=AJ{J+\D\) , (46) 

and therefore the breather frequencies in this case should 
lie above the linear band. 

The implementation of the AC limit can be achieved 
by setting J = and exciting one or several spins, so that 
they would start to precess around the hard axis with a 
frequency uj = 2|D|S'o, where Sq is the z projection of 
the spin. If the non-resonance condition lu ^ lOt^ — is 
satisfied, breather solutions can be continued. 

Breathers do not ff^^f^ in the continuum limit with 
easy-plane anisotropyllao. The reason is again that the 
corresponding linear band is gapless and unbounded, so 
that it covers the whole real axis. Correspondingly, there 
is no place for a frequency of a localised excitation on the 
real axis which does not resonate with the linear band. 
An essentially discrete modeled has been studied only in 
the case of a strong magnetic field directed along the 
hard axis. In this case the hard axis effectively becomes 
an easy axis and, as a result, the spins precess around 
the Z axis with a constant S^ component. In this case. 



a separation of variables ( 43 ) is possible which simplifies 
the treatment of the system. 



B. Breather solutions of the easy-plane ferromagnet 

As stated above, we do not apply an external mag- 
netic field, and therefore we do not change the ground 
state. As in the previous case of easy-axis, we compute 
breather periodic orbits from the AC limit using the gen- 
eralised Newton method. As a result we obtain solutions, 
schematically described in Fig. H (for one 'out-of-plane' 
spin) 




FIG. 6. Schematic representation of the discrete breather 
with two parallel 'out-of-plane' spins in the easy-plane ferro- 
magnet. 

For non-zero J previously non-excited spins start to 
precess with small amplitudes around the X-axis while 
the plane of precession of the 'out-of-plane' spin is no 
longer parallel to the easy plane, but is slightly tilted. 
Breathers with more than two precessing spins can be 
also created. 

Depending on its frequency, the breather width 
changes. When the frequency approaches the upper edge 
of the linear band, the breather becomes more delo- 
calised. However, this does not qualitatively influence 
its core structure, i.e., the effective precessing axis of the 
central spin is not continuously tilted towards the X axis 
upon lowering the breather frequency down to the linear 
band edge. The central spin dynamics can be viewed as 
a periodic (closed) orbit of a point confined to the unit 
sphere. Let the XY plane be the equatorial one. Then 
for large breather frequencies the point performs small 
circles around the north (or south) pole. Lowering the 
breather frequency does not change the fact that the loop 
still encircles the Z-axis. Thus the breather solution can 
not be deformed into a slightly perturbed and weakly lo- 
calised BEPW. This makes clear that the easy-plane fer- 
romagnet lattice supports breather solutions with a local 
magnetisation tilt which have no analog in a continuum 
theory. 

The situation is illustrated by Fig. M, 





n 

FIG. 7. Discrete breather profile J = 0.1, D = -1, 
^narrow = 1.1967 (diamonds) and uj^idc ~ 0.6649 (crosses). 
Inset shows the linear dispersion law and the location of 
breather frequencies. 

where the profiles of two breathers are represented: one 
corresponds to the frequency oJwide = 0.6649 , which is 
very close to the upper edge of the linear band w^ = 
0.6633 and the other one has the frequency tOnarrow = 
1.1967, which is far above the linear band (see the inset 
to Fig. |). 

The first solution is more delocalised, which can be 
seen from Fig. |7[ However, the central spin still pro- 
cesses in a way similar to the 'narrow' breather, i.e. it 
encircles the north pole on some lower latitude as com- 
pared to the narrow breather (see Fig. ||, two curves 
on the unit sphere correspond to two breather periodic 
orbits, discussed above). 




FIG. 8. Dynamics of the central spin with site number 
no = 15 in the chain, described in Fig. M. 

Hence, even when very close to the linear band, our 
breathers have a structure, which has no analogue in 
the continuum case. Moreover, we have investigated the 
dependence of the breather energy on the breather fre- 
quency (see Fig. O). We observe that there exists an en- 
ergy threshold since the breather energy attains a non- 
zero minimum when its frequency is still not equal to 
the edge of the linear spin wave spectrum. Note that 
for lattices of interacting scalar degrees of freedom dis- 
crete breathers have typically zero lower energy bounds 
in spatigiLdimension d — 1 and become nonzero only for 
d = 2, £1. The reason for the appearance of a lower 
nonzero bound in the present case is due to the already 
mentioned fact that the breather of the easy plane ferro- 
magnet system does not deform into a perturbed band 
edge magnon wave. Instead the central spin(s) is pre- 
cessing around the Z axis. This topological difference is 
the reason for the appearance of nonzero lower energy 
bounds. Such energy thresholds may be very impor- 
tant as they show up in contributions to thermodynamic 
quantities which depend exponentially on temperature. 
To eliminate possible size effects, we repeated calcula- 
tions for Fig. 1^ for a chain with A^ = 50 spins. The 
difference between the curves was negligibly small. 
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FIG. 9. Normalised energy E = H + JN/2 as a function 
of the detuning frequency f — uj — cj^ for a discrete breather 
with one out-of plane spin (curve 1) and with two out-of-plane 
spins (curve 2) with J = 0.1. The size of the system is A'' = 30 
spins. 

Energy thresholds can be estimated analytically in the 
limit of small exchange J. Ignoring displacements of 
all in-plane spins we obtain the threshold energy for 
a breather with M out-of-plane precessing spins, nor- 
malised to the ground state: 



E{uj) « 2M\D\S^ 



(47) 



where Sq = uj/2\D\ is Z-component of the precessing 
spin in the 

AC limit [see Eq. (0)]. The factor 2 comes from the 
fact that we should take into account the contribution of 
the breather tails. 

Equation ( |47| ) can be obtained using the following ar- 
gumentation. For large values of v in Fig. the main 
contribution to the breather energy comes from the AI 
'out-of-plane' precessing spins, because the tail ampli- 
tudes of the breather are small (see, for example. Fig. 
|7|). For 1/ — > the energy contribution from the tails is 
actually diverging. Thus the height of the minima of the 
curves in Fig. y can be estimated as two times the contri- 
bution coming from the central spins. Substitution into 
the above formula for the band edge frequency w^r yields 
E = E{ljj,) « 2MJ + 0{P). ] , For the case, considered 
in Fig. 0, for M — 1 (one precessing spin) the analytic 
result yields E « 0.2 while numerics give the value 0.33. 
In the case of two precessing spins (M = 2) the numerical 
result yields E « 0.44 while analytical estimate predicts 
the value 0.4. 

Increase of the frequency leads to a decrease of the 
precession radius of the central spin. In the AC limit, 
an upper bound for the breather frequency is defined by 
Lo = 2_D, which corresponds to the central (precessing) 
spin being parallel to the Z axis. This bound continues 
to exist when exchange is switched on. After reaching 
this frequency threshold, the breather becomes a station- 
ary (time-independent) solution. The existence of such 



a solution has been verified numerically by solving the 
time- independent Landau-Lifshitz equations. 



C. Stability of breather solutions and their 
asymptotic properties 

We have investigated the stability of our solutions with 
the help of the Floquet analysis (for details see Subsec. 
HID) and with direct Runge-Kutta simulations. It ap- 
pears that for small J's breathers with one precessing 
spin (see Fig. ||) are unstable, while a stable configura- 
tion which corresponds to two parallel precessing spins 
(see Fig. ^ is stable. Note that siroilar results have been 
obtained for the FPU-type latticcsc3. 

Stability tests also included the following nu- 
merical experiment. A periodic breather orbit 
{S'^{t)^°\si{t)^°\si{t)''°'^} is perturbed by deviating 
one of the central spins 5'^(0)(°) -^ S'fj(0)(°) +£ and sim- 
ulating the equations of motion (pi) . The error function 



A(t) = min , 
-e[o,T] ^ 



N 



n— 1 a— (a;,y,2:) 



(48) 



was calculated on each breather oscillation period T. In 
Fig. |l^, such a function (with e = —0.0025, in a chain 
consisting of A^ = 1000 spins with periodic boundary 
conditions) is shown for the breather solution with the 
two in-phase precessing spins (see Fig. ||). 
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FIG. 10. Time dependence of the effective error A for the 
breather solution with J — 0.1, D = —1 and oscillation period 
T = 4.4248. 

The error function is bounded during a significant pe- 
riod of time (more than 10000 breather oscillation peri- 
ods) and the breather structure remains preserved. Sim- 
ilar experiments have been performed for other types of 
breathers. They yield similar results. This demonstrates 
the stability of the discussed excitation. 

For a better understanding of the internal breather dy- 
namics, the data will be represented using the Fourier 
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expansion of the breather 
periodic orbit 



S^{t) = C^iuo- n) + ^ K(^; n) coskujt+ 

k=l 

+ B^{uj;n) sin kujt] , a = x,y,z. (49) 

We plot the space configuration of Cq and C„ = 
VAl + Bl 

The 'logarithmic' profile of such a solution is given in 
Fig. 11, We have plotted the space dependence of its 
Fourier harmonics (from the zeroth to the fifth one) for 
one particular stable solution. 

Let us analyse the behaviour and the exponential spa- 
tial decay of these harmonics. As can be seen from Fig. 
[ll| , the zeroth (static) component is present. According 
to this figure, the zeroth component decays exponentially 
in space. This seems to be surprising, since the corre- 
sponding zero frequency resonates with the bottom of 
the acoustic- type linear band [see Eqs. (|45|)-(|4q)]. 





FIG. 11. Spatial dependence of the Fourier components 
of the discrete breather of the type depicted in Fig. |8| for 
J = 0.18, D = —1, LO = 1.42. Numbers on the panels repre- 
sent the order of the harmonic k [see Eq. (US) for explanation] . 

To understand the results in Fig. n^ we linearise the 
equations of motion (^) around the ground state (^4|) in 
the breather tails. We obtain the following equations for 
5^ and S^ components {S^ is assumed to be equal 1 with 
higher than linear corrections): 



2\D 



- 2o,„ -f- o,,,_i 



-f5^_i) + 25,1 = 
= . 



(50) 



The numerical results suggest that the static S^ com- 
ponent is zero. This satisfies the second equation in (pfl) . 
The first equation in ( pO[ ) allows for an exponential de- 
cay of the static S^ component. Its decay can be charac- 
terised by the value Aq if Co(w; n) ~exp(— Ag|n|) , \n\ —> 
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oo , Ag > 0. The substitution of this ansatz into ( p5| ) 
yields 



(51) 



The spatial decay of all the other (non-zero) harmonics 
of the breather solution can be obtained from the disper- 
sion law (Ma) by substituting q = n — iX^ and solving this 
equation with respect to A with frequencies flk — kuj. As 
a result, we get 



K 


==ln 


L2W 

tj 



XI 



In 



c+ve 



1 



c 



VD^ 



^l \D\ 



J 



1 



1,2, 



(52) 



Since the Fourier components for S^ decay in space as 
S^ (except for the static one), we have omitted them. 
The spatial decay of the Fourier components of 5^ can 
be obtained using the normalisation condition (0), and 



therefore, for small deviations from the ground state (44), 
the following expansion is true: S^ = 1-3^/2-3^/2 + 

O ( 3n'^ ) • Substituting here the Fourier expansion for 

3^ and 3^ , we see that only the product of terms contain- 
ing harmonics kuj and {m ± k)uj of ^^'^will contribute to 
the decay of the mth harmonic of S"^. We have to choose 
the smallest exponent of all possible ones to obtain the 
leading order decay rate: 



/c— 0,l,...,m 



(53) 



As a result, the following relations have been obtained 
for the first five harmonics of the 3^ component: Aq = 

ZAi, Ai — Ai + A2, A2 — ^Aj^, A3 — A^ i- A2, A4 — 

Af -f A| . The comparison of these theoretical results with 
the values of A extracted from the numerical data is given 
in Table I. 



TABLE I. Numerically and analytically computed values 
of the decay exponents Xk of Fourier k = 0,1, .., 4 harmonics 
for different spin components. Parameter values are the same 
as in Fig. H. 



Order, k 


> 


X 

k 


K 






numerical 


analytical 


numerical 


analytical 





3.5914 


3.5903 


3.1846 


3.1856 


1 


4.9463 


4.8055 


1.8067 


1.7951 


2 


3.5859 


3.5903 


2.9979 


3.0103 


3 


4.8499 


4.8055 


3.5597 


3.5690 


4 


5.4393 


5.3641 


3.9195 


3.9309 
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As can be seen from this Table, the agreement between 
the numerical and analytical values of A decreases with 
the order of the Fourier components, which we think is 
due to the smallness of the higher order components. 



V. TWO-DIMENSIONAL LATTICE WITH 
EASY-PLANE ANISOTROPY 

Finally, we briefly consider a two-dimensional sys- 
tem, namely an easy-plane ferromagnet with nearest- 
neighbour exchange interactions. We have numerically 
simulated the Landau-Lifshitz equations for this system: 



qx — J qz ( cv 

'-'mn — '^y'-'mn \'-'m- 


-l,n + '-'m+l,n + '-'m,n-l + '^m.n+l) 


— JzSn \S,n-i 


n + "S^m+l.ri + ^m,n-l + "^^rn.n+lj 


~ '^D^mn^mn 




^mn = J^^mn \^m- 


-l,ri + ^m+l,n + ^m,n-l + ^rn,n+l) 


— Jxi>rnn \^m- 


-l,n + ^m+l,n + ^m,n-l + ^m,n+l) 


+ 2DS^„,„S:^^ 


(54) 


^mn — 'Jxt>rnn \^m- 


-l,n + ^m+l,n + ^m,n-l + ^m,n+l) 


~ '^V^mn \^m- 


-l,n + ^7n+l,n + '^m,n-l + "^m.n+lj 



using again a fourth-order Runge-Kutta scheme with var- 
ious initial spin configurations. The results of our simu- 
lations to some extent are similar to the one-dimensional 
problem. In Fig. |llj, we show the simplest possible con- 
figurations of breathers 
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FIG. 12. Schematical representation of some possible con- 
figurations of discrete breathers with four precessing spins. 



which involve four 'out-of-plane' precessing spins. No 
stable breathers with one precessing spin are possible, 
similarly to the one-dimensional model. Also, there are 
no stable breathers with two or three precessing spins 
(at least in the limit of small J). Among the three pos- 
sible stable configurations shown in Fig. n2, the first one 
[see panel (a)] corresponds to four spins precessing par- 
allel and in phase, similar to its one-dimensional coun- 
terpart. The second two configurations do not have ana- 
logues in the one-dimensional case, but have also similar 
parity properties. The cases shown in Fig. |l2|(b,c) rep- 
resent breathers with two spins precessing around the Z 
axis in the positive direction (marked by dots) and two 
spins precessing around the negative direction (marked 
by crosses). 

Simulations have been performed on a lattice of 150 x 
150 spins with J — 0.11 and D = —1. 



VI. SUMMARY AND DISCUSSIONS 

Summarising, we have considered discrete breathers in 
different (easy-axis and easy-plane) classical ferromag- 
netic spin lattices as essentially discrete objects. We 
have shown systematically how to implement the anti- 



13 



continuum limit for different types of magnetic lattices 
(different types of anisotropy). Depending on the type of 
anisotropy, discrete breathers have properties similar to 
breather solutions for other nonlinear lattices. In the case 
of an easy-axis anisotropy the breather solution appears 
in the gap of the linear magnon band as do the breathers 
of the Klein- Gordon- type models. In the easy-plane case, 
there is no gap in the linear band and the breather fre- 
quency lies above the band; these breathers resemble the 
breathers of the FPU-type chains (also known as the 
Sievers-TakenoQ modes). 

The concept of the anti-continuum limit helps us, first, 
to show rigorously the existence of discrete breathers, 
and, second, to compute the breather solutions numer- 
ically. The existence proof has been performed for the 
one-dimensional XYZ Heisenberg ferromagnetic chain 
with a strong exchange {Jx.y *C Jz)- The proof can 
be easily generalised to the presence of an easy-axis ion 
anisotropy and to larger lattice dimensions. The numer- 
ical 

continuation of the solutions from the AC limit has 
been done with the, help of a Newton iteration scheme. 
Note that so faJla only breathers with one nonzero 
Fourier component in time (Q) have been studied, due 
to the fact that it is much easier to treat them both nu- 
merically and analytically. The solutions we have stud- 
ied allow the infinite number of harmonics, as in XY Z 
model, for example. 

Why is it important to study discrete systems if the 
continuum approximation can give an analytical solu- 
tion? First of all, it is known that thfi^breathers are non- 
generic for most continuous modelstJ. Therefore many 
systems may be incorrectly referred to as those which 
do not possess breathers. We demonstrated this circum- 
stance for the easy axis ferromagnet, where the slightest 
exchange anisotropy in the hard plane leads to a loss of 
breathers in the continuum model, but not in the case of 
a spatial lattice. In addition we have obtained breather 
solutions for easy plane ferromagnets which have sim- 
ply no continuum analog. This is due to the fact that 
the spins in the center of the excitation precess around a 
tilted axis leading to a local tilt of the magnetisation. 

Finally, we would like to address some important unan- 
swered questions in this area. The first problem is how to 
treat quantum spin lattices (e.g., when the total spin is 
too small to treat the lattices classically) and what is the 
quantum analogue of the spin breather. Another impor- 
tant question is the breather's mobility. So fatyjthere is no 
rigorous existenp; proof for moving breatherscil, however, 
Lai and Sieverst£l have obtained some numerical results 
with highly mobile spin breathers. Since their results are 
concerned only with breathers with one Fourier compo- 
nent in time, it is still questionable whether breathers 
with an infinite number of harmonics can freely propa- 
gate along the lattice. 
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